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WILD  BEARINGS 

INTRODUCTION 

A  line  of  bearing  may  be  erroneously  attributed  to  a  particular  emitter. 
These  bearings  are  referred  to  as  wild  bearings. 

Several  algorithms  to  estimate  the  location  of  a  particular  emitter  and 
give  an  estimated  error  ellipse  (EEP)  utilizing  lines  of  bearing  (LOB's)  to 
the  emitter  have  been  proposed  in  the  literature.  The  possibility  that  some 
of  the  bearings  are  misidentified  (wild)  has  not  been  taken  into  account  in 
the  derivations  of  the  fix  algorithms  seen  to  date  by  MARC.  Accounting  for 
this  possibility  requires 

(1)  Estimation  of  the  probability  that  a  particular  bearing  is 
wild. 

(2)  Incorporation  into  a  fix  algorithm  of  the  estimate  of  the 
probability  that  a  particular  bearing  is  wild. 

It  is  the  purpose  of  this  report  to  show  how  this  possibility  may  be  taken 
into  account  by  incorporating  the  probability  of  a  bearing's  being 
mis-identified  (i.e.  being  a  wild  bearing)  into  fix  algorithms.  The  methods 
of  algorithm  modification  suggested  here  have  two  related  properties  that  seem 
particularly  desirable: 

(1)  Increasing  probability  of  misidentification  increases  the  size 
of  the  error  ellipses. 

(2)  Bearings  which  are  more  likely  to  be  correctly  identified  are 
weighted  more. 

The  body  of  this  report  is  divided  into  sections,  as  follows: 

Section  1  is  a  brief  discussion  of  considerations  that  affect  the 
probability  of  a  particular  bearing's  being  wild.  Section  I  also  contains  a 
reference  to  a  paper  discussing  similar  problems  in  the  context  of  radar 
tracking  algorithms. 

Section  11  illustrates  two  approaches  to  using  the  probability  of  a  wild 
bearing.  The  first  uses  an  analogy  with  the  ellipse  combination  methodology 
to  illustrate  the  qualitative  impact  that  should  be  expected  from  wild 
bearings.  A  method  of  modifying  likelihood  formulations  is  suggested  by  this 
analogy.  A  second  approach  is  also  shown  which  uses  a  weighted  average  of 
estimated  locations.  These  two  approaches  are  then  compared. 


4239  VIA  PADOVA,  CLAREMONT  CA  91711 


Sections  III  and  IV  show  examples  of  these  modifications  for  the 
sequential  and  perpendicular  fix  algorithms  respectively. 

Section  V  is  related  to  Section  I,  and  presents  some  of  the  mathematics 
of  determination  of  some  of  the  factors  affecting  the  probability  of  a 
particular  bearing's  being  wild. 

Finally  there  are  appendices  in  which  notation  is  defined  and  various 
technical  issues  (such  as  the  mathematical  details  of  converting  an  X-Y 
covariance  matrix  into  a  weight  on  line  of  bearings)  are  addressed. 

The  materials  presented  here  are  merely  meant  to  give  an  introduction  to 
the  problem  of  wild  bearings. 

I.  PROBABILITY  OF  A  PARTICULAR  BEARING  BEING  WILD 

Various  factors  might  be  considered  when  estimating  the  probability  of  a 
bearing's  being  wild.  Analysis  might  involve: 

(1)  Number  of  emitters  within  the  angular  sector  (confidence  region) 
centered  about  the  observed  line  of  bearing  (LOB).  This  number 
might  be  determinable  but  is  more  likely  to  be  an  estimate  based 
perhaps  on: 

(a)  Emitter  density  in  the  zone  of  Interest. 

(b)  Area  of  the  zone  of  interest  lying  in  the  angular  sector. 

(2)  Quality  (confidence  value)  of  the  emitter  identification.  There 
are  several  considerations  that  might  be  important: 

(a)  Quality  of  match  numbers  generated  in  the  emitter 
identification  algorithm. 

(b)  Comparative  likelihood  between  the  identified  emitter  and 
other  candidates. 

(c)  If  an  estimate  of  target  location  already  exists  then 
statistical  distance  between  the  observed  and  expected  lines 
of  bearing  is  a  useful  indicator  of  tdiether  the  new 
observation  merits  consideration. 

The  idea  of  mixing  of  various  likelihoods  and  probabilities  to  determine  the 
probability  of  correct  identification  has  already  been  used  in  the  context  of 
radar  tracking.  See  the  work  of  Bar-Shalom.  (For  example,  "Tracking  In  A 
Cluttered  Environment  With  Probabilistic  Data  Association"  Proc.  <ith  Symposium 
Nonliner  Estimation,  1973  pp.  13“22).  Because  (2c)  above  requires  the  least 
assumptions  about  the  emitter  identification  algorithm,  we  use  it  to 
illustrate  how  likelihood  can  be  used. 

Assume  that  considerations  such  as  those  in  (1 )  above  have  generated 
a  probability,  P^,  of  an  identified  bearing's  being  wild.  The  philosophy 

behind  (2c)  above  is  that  further  adjustments  are  necessary  because  a  bearing 
that  is  far  from  previous  observations  is  more  likely  to  be  wild  than  one  that 
is  close  to  previous  observations.  The  'usual*  procedure  is  to  assume  that 
the  probability  P^  is  modified  proportional  to  its  likelihood. 

(i)  Likelihood  is  based  on  a  statistical  distance  between  the  old 
emitter  location  estimate  and  the  new  bearing. 

(ii)  Statistical  distance  is  based  on  two  forms  of  uncertainty.  These 

are  the  uncertainty  in  the  estimated  location  of  the  emitter  and  the 
uncertainty  in  the  new  bearing.  The  details  depend  on  the  forms  in 
in  which  the  normal  distributions  used  by  the  fix  algorithm  are 
expressed.  Technical  details  are  discussed  in  section  V. 
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(iii)  To  compute  the  proportionality  it  is  first  necessary  to  compute 

the  average  likelihood  in  the  region  in  tdiich  LOBs  can  be  accepted 
for  this  emitter.  (This  is  often  done  by  integrating  over  this 
region  and  then  dividing  by  the  area  of  the  region.) 

(iv)  If  is  the  modified  probability  of  being  wild  (after  accounting 
for  fix  considerations)  then,  according  to  the  likelihood  approach 
P^:1-Pm  as  Py«(average  likelihood) ;(l-Py)*(observed  likelihood) 


The  details  of  this  analysis  are  completed  in  Section  V.  If  P^  can  be 
estimated  and  only  an  average  effect  of  misidentification  is  desired  then  P^ 
may  be  used  directly  saving  a  degree  of  computational  complexity. 

II.  USING  THE  PROBABILITY  OF  A  BEARING  BEING  WILD  IN  ESTIMATES  AND  EEP’s 


A.  Approach  1-  Modify  The  Theory  Behind  the  Various  Fix  Algorithms 


Fix  algorithms  examined  to  date  are  derived  by  minimization  of  a  quadratic 
form  related  to  the  log  of  the  likelihood.  (In  this  memo  that  quadratic  form 
will  also  be  referred  to  as  statistical  distance.)  The  modification  to  the 
quadratic  form  may  be  inferred  by  reference  to  likelihood  reasoning,  but  the 
same  analysis  is  somewhat  easier  to  state  in  terms  of  the  ellipse  combination 
algorithm.  This  reasoning  is  presented  below.  The  conclusion  is  that 
covariance  matrices  ^ould  be  divided  by  the  probability  that  the  data  are  not 
•wild.’ 

Compensating  for  a  probability  of  wild  ellipses  in  the  combination  algorithm 
is  straightforward.  Furthermore  the  result  suggests  a  way  of  thinking  of  wild 
observations  which  is  applicable  to  the  least  squares  formulations  used  to 
derive  the  standard  fix  techniques.  The  following  properties  of  the 
combination  algorithm  will  be  the  basis  for  the  analysis: 

(1)  The  combination  is  Independent  of  order. 

(2)  If  n  identical  ellipses  are  combined  then  the  resultant  ellipse 
has  a  covariance  matrx  •  1/n  times  the  original. 


The 


Now  imagine  that  a  fraction  P  of  the  ellipses  are  expected  to  be  wild, 
correct  result  will  be  obtained  if  the  resultant  covariance  matrix  is 
multiplied  by  1/(1-P).  Since  combination  is  independent  of  order  the  same 
result  would  apply  regardless  of  the  number  of  individual  covariance  matrices 
mixed  in. 

Thus  given  a  theory  which  is  based  on  minimizing  a  least  squares  function 
of  emitter  location  (X,Y)  (equivalent  to  maximizing  the  normal  likelihood 
function)  such  as  the  minimization  of  squared  angular  error 


j(x,Y)-^jy  CVe,^(x.Y)fp2c^en-e^(x.Y)] 
there  is  a  simple  modification  to  account  for  an  estimate  of  the  probability 


of  being  a  wild  bearing.  If  is  the  probability  of  the  kth  bearing  being 


wild  then  the  modified  version  is 

jp(x,Y)-^5y  (i-p,j)Ce^-e^(x.Y)3^p^[e^-e,^(x,Y)] 


Derivation  of  location  estimates  and  EEP's  (error  ellipses)  from  this 
form  is  straight  forward  in  that  the  old  algorithms  are  used  with  the  small 
modification  that  formulas  involving  observation  covariances  are  divided  by 
the  probability  of  that  observstion*s  being  true. 


B.  Approach  2-  Use  a  weighted  average  of  the  original  estimate  and  the  update 
that  would  apply  if  the  new  bearing  were  known  to  be  true. 


The  weight  used  is  the  probability  that  the  bearing  is  true.  If 

S^«the  estimate  and  covariance  matrix  without  the  new  bearing 

S2*the  estimate  and  covariance  matrix  with  the  new  bearing 

(x.y).  S  >the  combined  estimate  and  covariance 


Then  for  i-pz, ♦(1-p)22  it  is  shown  in  Appendix  3  that 

2  ^  2 

x.y.  X,  x-y- 

S-PCS,^(  ^  ^  l)>{ 


Note:  This  approach  is  analogous  to  one  found  in  the  Bar-Shalom  paper 
previously  referenced. 


C.  Comparison  of  approaches  1  and  2. 

Approach  2  only  works  well  in  its  current  form  for  recursive  algorithms. 
See  for  example  the  approach  2  comments  for  the  perpendicular  method  (which  is 
not  recursive)  described  in  section  IV.  There  are  theoretical  considerations 
which  suggest  that  when  approach  2  applies  it  is  superior.  These 
considerations  are  discussed  below. 

Approach  1  only  accounts  for  one  of  two  ellipse  size  considerations, 
however.  If  p«.1  for  each  of  10  LOBs  added  then  approach  1  is  similar  to 
adding  9  LOBs  with  p-0.  Approach  1  does  not  make  any  adjustment  for  the 
Increase  in  uncertainty  owing  to  inclusion  of  a  wild  bearing.  However, 
approach  2  makes  such  an  adjustment  and  is  generally  more  robust  (i.e.  it 
accounts  for  some  error  types  not  built  into  the  model).  The  covariance 
matrix  generated  by  approach  2  has  more  ways  to  slow  EEP  shrinkage  if 
innovations  are  bigger  than  expected.  Simulation  has  not  yet  been  done  to 
check  the  implications  of  the  theoretical  considerations  concepts  discussed 
here. 

Further  investigation  may  reveal  some  method  of  combining  these  two 
approaches  to  get  the  advantages  of  both. 


III.  WILD  BEARINGS  AND  SEQUENTIAL  ESTIMATION 


Notation  is  defined  for  what  follows  in  the  notation  appendix.  It  will  be 
assumed  that  the  probability  of  a  wild  bearing,  is  known. 

The  remaining  portions  of  the  algorithm  for  sequential  estimation  will  be 
discussed  using  both  of  the  approaches  of  Section  11. 

A.  Approach  1  -  The  Likelihood  Modification  Method 
The  modified  objective  function  using  approach  1  is 

The  update  for  the  location  estimate  becomes 


The  natural  choice  for  the  covariance  of  the  estimate  is 


5- 


B.  Approach  2  -  Weighted  Location  Estimate 
^®t  ^‘*k*^k^*^*Bk-1*^0kk-1*ek^0kk-1  ^  ^*^0kk-1*0k^®k”®k^®k-1 


*Bk*®k^*Bk-1*^ 


‘k-1 

ViVl 


ViVi 


VI  "0kk’^0k  0kk' 


-InT  ^-1 


K. 


k-1 


Vk. 

“k^  -k 

Note  the  similarity  of  the  location  estimate  to  that  used  in  method  1. 

IV.  WILD  BEARINGS  AND  THE  PERPENDICULAR  METHOD 

Let  denote  the  location  of  the  detector  of  the  kth  LOB. 

Let  0^-the  angle  measured  clockwise  from  true  north  to  the  LOB. 

Let  M  '*^"®k®®»®kj 

-sln§  cosS.  sin^S, 

k  k  k 

Let  d^ ^ ^ ^ ^ **“k ^ ^ ^ ^sln^§l^-2( x-u^ ) ( y-Vj^ ) sin§|^cos0|^ 

■the  square  of  the  distance  from  (x,y)  to  the  kth  LOB 
(measured  along  a  perpendicular  to  the  LOB) 

It  will  be  assumed  that  the  probability,  P|^,  of  the  kth  LOB's  being 

wild  is  known.  Location  estimation  will  be  examined  using  the  two  types  of 
modification  discussed  in  Section  II. 

A.  Approach  I  -  The  Likelihood  Modification  Method 

The  function  to  be  minimized  without  wild  bearings  is  a  variation  of 

d®(x.y) 

The  variations  are  modified  weightings  in  this  sum.  The  wild  bearing 
modification  is  of  this  type.  The  modified  likelihood  is 

k£? 

The  resultant  point  estimate  is  at 

(x,y)^-[j^ry(1-Pj^)Mj^]  k5l^”k^“k’^k^^^ 

B.  Approach  IX  -  Weighted  Location  Estimate 

Ipdt  denote  the  Nth  location  estimate.  Then 

(x3.y3)^-P3(Xj,,^,yjj.l)^^(1-Pjj)C^5jM^]  ’ 

The  problem  with  this  approach  is  most  easily  seen  in  the  ease  where  the 
probability  of  a  wild  bearing  is  sufficiently  small  that  the  probability  of 
two  wild  bearings  is  negligible  and  can  be  dropped.  After  dropping  the 
negligible  terms,  are  lost*  To  resolve  this  problem  it  is 

neccessary  to  track  at  least  N  estimates  (and  ideally  2®)  depending  on  which 

estimates  are  wild. 
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V.  WEIGHT  (PROBABILITY)  OF  A  PARTICULAR  BEARING  BEING  WILD 

The  generation  of  a  statistical  distance  as  suggested  in  Section  I 
consideration  (2c)  above  requires  mixing  the  distributions  of  two  types  of 
error.  These  two  types  of  error  are: 

(1)  error  in  one's  previous  estimate  of  the  location 

(2)  error  in  observed  angle 

Because  these  errors  are  described  in  different  coordinate  systems  it  is  not 
easy  to  combine  them.  In  what  follows  we  make  the  same  simplifying  assumption 
that  was  made  in  deriving  the  perpendicular  method.  That  assumption  is  that 
perpendicular  distance  from  the  target  to  the  line  of  bearing  (LOB)  can  be 
used  instead  of  angular  error.  (For  a  partial  solution  of  this  problem 
without  making  the  perpendicular  distance  approximation,  see  the  Mathematical 
Appendix. ) 

It  is  necessary  to  decide  if  procedures  ensure  that  if 

instead  other  assumption  is  appropriate.  In  fact  it 

might  be  better  to  assume  until  it  is  certain  that  a  target  of  some 

sort  has  been  identified.  After  that  point  has  been  reached  we  assume  that 
two  proportionalities  affect  the  weighting: 

(likelihood  observed)  :  (likelihood  of  an  average  wild  bearing)  (1) 

1-Pw  :  (2) 

Evaluated  and  combined  this  yields  that  weights  are  proportioned  according  to 

where  N(a:b,c)*the  value  of  the  normal  density  with  mean  b,  variance  c, 
evaluated  at  a.  The  acceptance  value  referenced  would  typically  be 
independent  of  index  i. 

From  this  proportionality  relationship  and  since  the  sum  of  the  two 
probabilities  is  one,  we  have 

Pi-(»il/Lj)«P^/[N(0i-e(Bi.,);O,t|»ci)*(1-Pw)*(Pi/Lj)«P^^] 
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APPENDIX  1  *  RECTANGULAR  TO  ANGULAR  COVARIANCE  CONVERSION 

The  likelihood  which  must  be  converted  from  rectangular  to  angular  coordinates 
is  given  by 

L(x.y)-(1/(2ii|S|  )exp(-(x-XQ,y-yQ)S"^  (x-Xp.y-yQ)'^)  where 

•previous  location 
estimate 

In  this  appendix  a  density  L(0)  is  computed  from  this  density.  This  is  done 
in  section  A  below.  Since  the  result  is  not  normal  an  approximation  is 
computed  in  section  B  using  some  of  the  same  approximations  that  were  used  in 
derivation  of  the  perpendicular  method. 

A.  L(G)  Exactly. 

To  express  the  likelihood  function  in  polar  coordinates  first  define 

•1  T 

a-(sineiCose)S  (sine.cosG) 

b-(Xo,yo)S“^ (sine.cose)^ 

The  likelihood  function  in  polar  coordinates  is: 

f(r,0)r  dr  d0-(1/(2ir|S|  ))exp(-(ir  -2br+c))r  dr  d0 

Elimination  of  r  to  get  a  likelihood  depending  only  on  0  requires  integrating 
this  expression  with  respect  to  r.  First  complete  the  square  using  the 
following  terms. 

Let  R«b/a  D«c-aR^ 

g(0)-(1 /(2ir|S| ))  /•  exp(-(ar^-2br+c))r  dr 

•  (1/(2ir|S|))  /q  exp(-a(r-R)^-D)r  dr 

-(exp(-D)/(2v|Sj ))  /"jj  expC-aC*— R)^)((r-R)+R)  d(r-R) 

There  is  usually  no  probability  of  r<0  so  it  should  be  a  reasonable 
approximation  to  use 

(exp(-D)/(2x|Sj ))  /\  exp(-a(r-R)^)({r-R)+R)  d(r-R) 

Odd  functions  integrate  to  0  therefore  the  above  is 
-{exp(-D)/(2ii|S| ))  /"^  exp(-a(r-R)^)R  d(r-R) 


•  (exp(-D)/(2ii|S| ))  (i/iTz/a)  R 

It  is  necessary  to  multiply  this  density  by  a  normal  and  integrate  from  0 
to  It.  This  is  simplified  somevdiat  by  evaluation  of  d  which  turns  out  to  be 

CxQCOse-yQSine]^/(a|S| ) 

It  is  also  simplified  by  assuming  that  a  is  a  constant  (if  the  ellipse  were  a 
cicle  it  would  be.)  Recall  however  that  R>b/a.  Thus  %fe  reduce  to 


-b( exp( -C  XQCose'yQSine]^/(a I S | ) ) /( 2a | S | /wa ) 


B.  A  Normal  Approximation  for  L(e) 

In  this  approximation  the  first  step  is  to  evaluate  at  the  closest  point 
rather  than  integrate  over  all  points.  First  note  the  closest  point 
(x-XQ,y-yQ)"±(distance  from  LOB  to  location  estimate)(-cose(Bj)),sin0(B^)) 

■±r^sin(e-0(B^))(-cos(e(Bj)),sin(0(Bj))) 

-±rj(0-0(Bj))(-cos(0(B^)),sin(0(Bj))) 

This  form  of  the  approximation  is  just  what  is  needed  to  rewrite 
(x-Xq)s“^  (x-Xq)^  as  (©-0(Bj))i()”^  (0-0(B^)) 

In  particular  41  ^•r^(-cos(0(Bj)),sin(0{Bj)))S  ^ (-cos(0(B^) ,sin(0(B^ ) ) ) 


APPENDIX  2  -  NOTATIONAL  APPENDIX 


Let  P^-the  probability  that  an  average  bearing  is  wild.  (Given) 

Pj"the  'probability*  that  the  ith  observed  bearing  is  wild.  (Goal) 

©.•observed  angle  for  ith  LOB  measured  clockwise  from  North. 
(Observed) 

B^-estimate  of  target  location  after  i  bearings.  (Recursive) 

Bj  is  a  vector  (x^.y^). 

6^(Bj)*angle  from  detector  of  the  ith  LOB  to  B^.  (Computed) 

r^-dlstance  from  to  the  detector  of  LOB  1. 

♦Qj^«covarlance  matrix  for  the  observation  0^. 

-estimated  covariance  matrix  for  the  estimate  B. .  (Recursive) 

,  -  -1  -  .  T  ii 

♦ci“^f  CC"Cos(ei),sin(gi)]<)Bj_.,["COs(ei),sin(ei)]M 

See  the  Appendix  for  a  discussion  of  the  problems  in 
determining  ♦ci*  This  is  an  approximation. 

L  -Length  of  the  acceptance  interval  (in  0)  determined  by  . 

X  Lk  X 

Uj-probabllity  of  being  in  the  acceptance  Interval  (of  length  L^). 

j-[-eoa(0j  (Bj ) ) ,  sin(0j  (Bj ) )  ] 


APPENDIX  3  '  WEIGHTED  AVERAGE  COVARIANCE  DERIVATION 


Given  that: 

a)  An  estimate  z^-(x^,y.|)  would  be  used  if  the  latest  bearing  was  known 
to  be  wild. 

b)  An  estimate  Z2-(x2iy2)  would  be  used  if  the  latest  bearing  was  known 
not  to  be  wild. 


c) 


d) 

e) 


The  covariance  for  z^  would  be  (if  it  applied)  and  the  covariance 
for  Z2  would  be  if  it  applied.  Furthermore, 


The  probability  that  the  bearing  is  wild  is  p.  For  notational 
convenience  also  define  p^*p  and  P2~1~P> 

The  combined  estimate  z«p^z^^P2Z2“Pz^*(l-p)Z2. 


The  combined  estimate  z  has  covariance  § 
-  f? 


p^f\/2^(2-z)(2~2)^S(z;Z^»S^)dxcly 


■i£f  PiCC(,. 


(x-x)*^  (x-xXy-y) 


(x-x)(y-y)  (y-y)^ 

Thus  the  i  entry  of  S  is 

^“iSl  N(z;rj,Sj)dxdy 

1 

-  Pj  [aj>2x^-x^-2xj 

•  iJT  Pi 

-  p^  Ca^  ♦(x-x^  )^>P2Ca2*(x-X2)^D 


)  N(z;Zj,Sj)dxdy 


ill  Pj/r.J’r,C(Xi“x)^*2xXj-Xj'2xx+(x)^]N(z;z^,Sj)dxdy 
^  2  „2 


-  p^Ca^>x|]>pp2Ca2^x|]-(x)^ 


Similar  expressions  are  obtained  for  the  other  entries  of  S.  The  combined 


as  X  y 

S-p^CS,^(  ^  ^ 


*2^2 


)>( 


(x)2 


X  y 


x,y, 


x-y 


2'2 


X  y 


(y) 


In  some  cases  it  might  be  preferable  to  use  the  equivalent  form 


(x^-xg)^  (x^-XgXy^-yj), 


S-pS,  +  (1-p)S2^p{1-p)(  ’  ^  ^ 

(x^-x2)(y,-y2)  (y,“y2)^ 


